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Abstract 

We consider the spectral and dynamical properties of quantum systems of n 
particles on the lattice Z**, of arbitrary dimension, with a Hamiltonian which in 
addition to the kinetic term includes a random potential with iid values at the lat- 
tice sites and a finite-range interaction. Two basic parameters of the model are 
the strength of the disorder and the strength of the interparticle interaction. It 
is established here that for all n there are regimes of high disorder, and/or weak 
enough interactions, for which the system exhibits spectral and dynamical local- 
ization. The localization is expressed through bounds on the transition amplitudes, 
which are uniform in time and decay exponentially in the Hausdorff distance in 
the configuration space. The results are derived through the analysis of fractional 
moments of the n-particle Green function, and related bounds on the eigenfunction 
correlators. 
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1 Introduction 



1.1 On localization in the presence of interactions 

In the context of non-interacting particle systems, or equivalently one-particle theory, 
Anderson localization is a well studied phenomenon, which for various regimes of the 
parameter space can be established even at the level of rigorous mathematical analy- 
sis (e.g. ICL90, PF92 StOl , ASh-01 , GK06, Ki08b | and references therein). The picture 
is far less complete when it comes to systems of interacting particles subject to a ran- 
dom external potential, which generally may be expected to produce localization. 

Particularly perplexing is the situation where there are n fermions in a region of 
volume |A|, with |A| — > oo and n/\A\ p > 0. It was proposed, through reasoning 
presented in PBAA061, that if the interactions are week and the mean particle separa- 
tion is significantly below the localization length of the non-interacting system, then 
the interaction would not affect by much the dynamical properties of the system. 

In particular, such reasoning has lead to the suggestion that if the system is started in 
a configuration for which the density of particles in one part of the region is higher that 
in another then the uneven situation will persists indefinitely, assuming the Hamiltonian 
is time independent. While that would be in violation of the equipartition principle, 
it would be in line with the dynamical behavior of the the non-interacting system in 
the regime of complete Anderson localization for the one-particle Hamiltonian ( IIAi94l 
IAS-hOIIIGKOTI ). 

Rigorous methods are still far from allowing one to decide whether complete lo- 
calization will persist in the presence of interactions, as claimed in IIBAA06II . Fur- 
thermore, the analysis of even a fixed number of particles with short range interactions, 
and |A| oo, has presented difficulties. An important step was recently made by Chu- 
laevsky and Suhov IICS07IICS 08b | who proved the existence of spectral localization for 
systems of two interacting particles which are subject to highly disordered external po- 
tential. The authors expect that the analysis of the n = 2 case, which is based on the 
multiscale approach of IIFS83I|DK89I , could be extended to any finite n. In this work 
we approach the question using somewhat different tools, and address also the issue of 
dynamical localization. We establish the existence of localization regimes for any fi- 
nite N, with decay rates which are uniform in the volume. Curiously, as is indicated in 
the figure below, the bounds which are established here carry a qualitatively somewhat 
stronger implication for n = 2 than for higher values of n. 

1.2 Statement of the main result 

Our goal here is to present a basic proof of localization for an arbitrary number of 
particles moving on a lattice of arbitrary dimension, which for convenience is taken to 
be Z"*, in regimes of high disorder or sufficiently weak interactions. 

There are a number of ways to formulate a quantum system of particles on a lat- 
tice, which here is taken to be Z''. We find the following convenient, but the method 
discussed here can be also be adapted to other formulations. 

The Hilbert space of n particles on Z'^ is the direct product H^") = e^{Z'^y\ We 
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take the Hamiltonian to be an operator of the form: 

n 

H^'^\uj) = ^[-Aj + AF(2;,;c^)] + Zi(x; a) , (1.1) 

with: A the discrete Laplacian (second difference operator) in l/, V{x^ lo) a random 
potential (described below), and U a finite range interaction which is given in terms of 
functions of the occupation numbers 

p 

U{-k;cx) = J]"'^ Y1 Ua{Ma{-^)) (1.2) 

where 

AA4(x) - {iV,(x)}„g^, (1.3) 

with Nu{:x.) — ^xj,u the number of particles the configuration x e (Z'^)" has at 

u g Z'*. It is to be understood that Ua{x) = unless Iltieyi ^n(x) ^ 0. 

The family of Hamiltonians is parametrized by A G M^, which controls the strength 
of the disorder, and a. :— {a2, ■ ■ ■ , Oip) G M^^^ which is the strength of the interaction. 
Obviously, the value of ak is of relevance for H^^\uj) only for n > k. It will be 
assumed throughout the paper that: 

Al The random potential is given in terms of a collection of iid random variables, 
with 

E [exp(t 1 1/(0)1) ]< oo foralHeK, (1.4) 
whose probability distribution is of bounded density, i.e., 

F(yix) e d^) = with g e L°° , (1.5) 

satisfying: 

g{v) < K Q{v~u)du, forallMGM, (1.6) 

J\v'\<Eo 

at some Eg < oo and i^T < oo. 

A2 The interaction terms are bounded, with |[/^(n)| < 1 for all A C Z"^ and all 
n G (Z'*)l"*l, and translational invariant, i.e., Ua = Ua' if A' is a translate of A. 

The above assumptions could be relaxed. In particular, translation invariance can 
be replaced by suitable translation invariant bounds, and, as in the case of one particle 
localization theory, the absolute continuity of the measure and ( |1.6| l can be replaced by 
a local power-law concentration bound such as the following condition: 

V{V{x) e[v + e,v-e]) < e^Kr V{V{x) e[v + Eo,v- Eo]) . (1.7) 
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for some t G [0, 1) and all < £ < 1. Under such reduced assumption, for which the 
case T — 1 corresponds to ( |1.6| l, the fractional moment bounds presented below would 
be limited to < s < t (and minor adjustment will be required in the argument, cf. 
BAS+Olll ). but that would not adversely affect the main results. 

Our main result is naturally stated in terms of the eigenfunction correlator which is 
introduced in Section |4] However, the statement can also be presented as follows. 

Theorem 1.1. Under the assumptions Al - A2 , for each p and n G N there is an open 
set in the parameter space, v'i^^ C M-|_ x M^^^ such that: 

1. For all (A, a) G r^f-* and up to n particles, i.e., k G {1, . . . , n}, each operator 
H^^^iij) has almost surely only pure point spectrum, with the corresponding 
eigenfuntions being exponentially localized in the sense o/dist7^(x, xq), as is 
explained below. 

2. Furthermore, for all (A, a) G T^f' , all k G {1, . . . , n}, and all's., y G {'L'^)^: 



E 



sup 

/lloo<l 



< Ae-'^'"*"^^'^)/^ , (1.8) 



where 



dist-M(x,y) := max max dist({a;i}, Y), max dist({j/i}, X) > , (1.9) 

\l<i<k l<i<k J 

with Sx"^ ['^y^''] the k-particle position eigenstates of corresponding to x [y], and 
the constants A and ^ depending on n, p but not on k (< n). 

3. The localization region T^"' includes regimes of strong disorder and of weak 
interactions, i.e., 

(a) for each a G W^~^ there is X{a) such that T^'' D (A(a), oo) x {a}. 

(b) for any A G t\^\ i.e. one for which the one-particle Hamiltonian exhibits 
complete localization, there are a{X)j > , j = {1, such that 
includes all {X, a') for which \a'\ < |a(A)| componentwise. 

The bound ( |1.8| l, applied to f{H) = e^'*^, implies dynamical localization, and 
through that also the spectral assertion which is made in Theorem |1.1| The latter is 
explained more explicitly in Appendix [B| 

One may note that the distance between configuration, dist-H(x, y), which appears 
above corresponds to the Hausdorff distance between the sets X = UiLii-^*} 
Y = Ur=i{2/i}' ^^^^ subsets of Z'' with its Euclidean metric. The exponential 
bound presented above deserves a number of further comments. 
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1.3 Remarks on the rate of exponential decay 



For systems of non-interacting particles, i.e., the case cx — 0, the one-particle localiza- 
tion theory allows to prove that in regimes of localization a bound like ( |1.8[ ) holds with 
a stronger decay rate. For the stronger bound (which can be established for regimes of 
strong enough disorder, or extremal energies, and in one dimension the full range of 
energies) the relevant distance is not distH(x, y) but: 

n 

dist(x,y) := 5]|xj-yj|. (1.10) 

One could ask whether exponential decay in dist(x, y) persists also in the presence 
of interactions. Upon reflection, the general answer to this should be negative: The n 
particle Hamiltonian clearly commutes with elements of the permutation group 
Sn- In the non-interacting case (i.e., a ~ 0) its spectrum is degenerate, H^'"\lj) 
being the sum of n commuting unitarily equivalent operators, each affecting only one 
particle. However, since generically the interactions couple the different permutation- 
related degenerate eigenstates of the non-interacting systems, it is natural to expect 
that for a ^ 0, the operator H'-'"\uj) will have no eigenstates in which the probability 
amplitude is essentially restricted to the vicinity of a particular n particle configuration. 
Instead, localization may still be manifested in the existence of eigenstates which decay 
in the sense of the symmetrized configuration distance 

n 

dists(x,y) := min V |xj - y^j | , (1.11) 

with Sn is the permutation group of the n elements {1, n}. 

The dynamical version of this eigenfunction picture is that for very large t a state of 
the form e~**^' ^ Sx^\ which has evolved from the initially localized state at x, would 
have non-negligible amplitude not only in the vicinity of x but also in the vicinity of 
the permuted configurations ttx. 

The above considerations are of course superfluous in case one is interested only in 
the fully symmetric or antisymmetric sector, where the initial states cannot be localized 
in the stronger sense and where only decay rates which are symmetric under permu- 
tations are of relevance. However, the decay rate exp(— dist-H(x, y)/^) is still qual- 
itatively weaker than exp(— dists(x, y)/^'). In particular, for configurations which 
include some tight subclusters with multiple occupancy our bounds do not rule out the 
possibility, which we do not expect to be realized, of the excess charges being able 
to hop freely between the different subclusters, as depicted in Fig. 1. Nevertheless, 
the bounds allow to conclude the main features of localization. Some further explicit 
comparison between different distances are given in Appendix [A| 

1.4 Comments on the proof 

It would be natural to ask why is there a need for a separate proof of localization 
for the n-particle system, since the configurations of the system can also be viewed 
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Figure 1: A schematic depiction of our localization bounds: Starting from the config- 
uration depicted on the left, at any later time except for events of very small proba- 
bility the collection of locations which are near an occupied site does not change by 
much. However, the Hausdorff metric bounds allow for the possibility that if the initial 
configuration had two or more particle in sufficient proximity [within the localization 
distance] then such 'excess' may transfer among the occupied regions. 



as describing a single particle with nc?-dimensional position vector, x = {xi, x„). 
Regarded from this perspective, the Hamiltonian may at first appear to have the 
usual structure for which localization is well understood, consisting of the usual kinetic 
term, a potential function C/(x), and a random potential J^j ^i^j'-:'^)- Th^ answer 
is that the values which the random potential V{x; ui) assumes at different positions 
in the nd-dimensional space are not independent. Instead, they are correlated over 
arbitrary distances, and the number of its independent degrees of freedom {L'^), for the 
systems in a box of linear size L, scales as only a fractional power of the number of 
its configurations (L""*). From this perspective, the randomness is much more limited 
than what is found in the well understood one-particle situation. 

The proof of Theorem |1.1| is organized as induction on n, and is guided by the 
following picture: once localization is established for less than n-particles, one may 
expect that throughout most of the volume the time evolution of n particles does not 
disperse, except possibly when the particles are all close to each other and move as 
some n-particle cloud. This hypothetical mode, with the n particles particles forming 
a quasi-particle, is ruled out using one-particle techniques, which are modified to show 
that such a possibility does not occur at weak enough interactions. 

Technically, our proof makes use of the Green function fractional-moment tech- 
niques, and in particular the finite-volume criteria of fAS+OT^. However, in addition 
to adapting a number of "off the shelf" one-particle arguments we need to show that 
exponential decay of the fractional moment of the Green function for lower numbers 
of particles implies: /. exponential decay for systems composed of non-interacting 
subsystems, and //. for the interactive system - exponential decay in a distance defined 
relative to the set of clustered n-particle configurations. These terms are explained 
more explicitly in the following sections. 

In Sections|2]-[5]we present a number of relations which are utilized in the derivation 
of Theorem | 1 . if These are then strung to a proof in Section|6] 



2 Finiteness of the Green function's fractional moments 



The proof of localization proceeds through exponentially decaying estimates for the 
Green function Gq{x, y; z) of finite volume versions of H^^'^\ijj) evaluated at energies 
z within the spectrum of the infinite volume operator Our first step is to establish that 
for s < 1 each |Gsi(x,y; z)Y is of finite conditional expectation value, regardless of 
z e C, when averaged over one or two potential variables - provided each of the con- 
figurations (x and y) includes at least one averaged site. The basic strategy is familiar 
from the theory of one-particle localization. However the proofs need to be revisited 
here since we are now dealing with random potentials whose values for different con- 
figurations are no longer independent, and in certain ways are highly correlated. 

It may be noted that the celebrated Wegner estimate is not being explicitly used 
in the Fractional Moment Analysis, its role being taken by the finiteness of the Green 
function's fractional moments. However, in view of the Wegner's estimate's intrinsic 
interest and conceptual appeal and we comment on it in Appendix [P] 

We shall use the following notation: The n-particle Green function associated with 
some region Q. Ql/^ and z E C+, is 

Go(x,y;z)EEG[,")(x,y;z):-^4"),(ij(f)-z)"'4")^ . (2.1) 

where iJ^"'' (w) is the (finite-volume) operator obtained by resticting to the Hilbert 

space := i?^(ri)" (with the default choice of Dirichlet boundary conditions). The 

vectors 4"'' , <5y"'' G ■* correspond to localized states, i.e. (4"'' , ^) = ^/'(x), and are 
parametrized by configurations x = (xi, . . . , a;„) of n-particles. When clear from the 
context we shall drop the superscript (n) at our convenience. 

The set of configurations with all particles in fl is denoted by C^"^ (fi) :— fi". Also: 



For a given set S* C fi, we denote by C''"'(ri; S) the set of n-particle configura- 
tions which have at least one particle in S. In case S = {x} the set will also be 
denoted asC^"^^;^^)- 

We denote by Cr^\fl) := {x e fi" | diam(x) < r} the set of configurations 
with diameter less or equal to r, the diameter of a configuration being defined as 

diam(x) maxj.fc jaij— Xfcj. 

Theorem 2.1. For any s € (0, 1) there exists Cs < oo such that for any ft C Z'', 
any two (not necessarily distinct) sites Ui,U2 G f2 and any pair of configurations 
X g C^"-'(f2;ui) andy € C*-"-'(ri; U2), the following bound holds 



E 



(|G[r^(x,y;z)[ | < G. (2.2) 



for all z £ C+ and A 7^ 0, with ^ — 2 in case ui ^ U2, and ^ — I otherwise. 
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For the proof, let us note that in its dependence on the single-site random variables 
V{x) the Hamiltonian is of the form 

Hn{io) = A+XY,V{u;Lo)N^, (2.3) 

where iV„ is the number operator, {Nuip) (x) := ^xk,u V-'(x), which counts the 

number of particles on the site m G Z'^. 

In analyzing averages over the potential variables we shall employ the follow- 
ing double sampling bound, which is the dual form of the regularity assumption Al, 
Eq. ( |1.6| l. Under that assumption, for any non-negative function h of one of the single 
potential parameters V = V{u), for some u e Z'': 

h{V) p{V) dV < KEo f f h{V + V) p{V) dV . (2.4) 

JrJ\V'\<Eo ^0 

Proof of Theorem \2.1\ We shall estimate the conditional expectation in \2.2\ with the 
help of the double sampling bound ( |2.4| l , applied to the pair of variables V{uj) (or 
a single one in case they coincide). According to that, it suffices to estimate just the 
integral over the variable(s) V'{uj) of the fractional moment of 

G^(x,y;z) := (<5,,(i/+ V {w)N^ ~ z)'^ 5y) , (2.5) 

At this point, a useful tool is the following weak L^-estimate which forms a straight- 
forward extension of ||AE-h06I Prop. 3.2]: For any pair of normalized vectors t/j in 
some Hilbert space, any pair of self adjoint operators with M > 0, and a maximally 
dissipative operator K: 

1 C 
> d^dTj< J, (2.6) 



1 



for all t > with some (universal) constant C < oo, where 1 [. . . ] denotes the indicator 
function. A similar bound holds for the one-variable version of ( |2.6| ). 

Applying ( |2.6| l to the expression in ( |2.5| l, and noting that 6^ and 6y are eigenvectors 
of Nu-^ and N^^ with eigenvalues greater or equal to one, one gets: 

W{t) E^^\{{V'{ui),V'{u2)) e[-Eo,Eof ■.\G'ni^,y;z)\>t}\ 

where W{t) is introduced just for the next formula. To estimate the corresponding 
integral of the kernel's fractional moment, one may use the Stieltjes integral expression: 

|Gn(x,y;z)| — = / W{t) d{t ) 

l-Eo.E„]^ ^0 -C'O Jo 

- ii-s)i\\\Eoy ^^-^^ 

A similar bound holds for in case ui — U2 for the average over a single variable. The 
bound (|2.2[) is implied now through a simple application of \2A\. □ 



9 



3 Localization domains in the parameter space 



The following notions are useful in describing localization bounds which persist when 
the strength of the disorder is driven up, and also to present the localization regimes 
which are discussed in this work. For their formulation we denote by Ai G M+ the 
critical coupling above which the one-particle Hamiltonian H^^^ exhibits uniform 1- 
particle localization in the sense of Definition 3.3 below. Its existence was established 
in flAM93L 

Definition 3.1. . A robust domain in the parameter space is a non-empty open set 
r C M+ X RP-i for some peN such that: 

1. if (A, a.) e r, then for all A' > A also (A', a) eT 

2. for every a e M^^^ there exist A(q;) G R+ such that (A(q;), a) G F, 

3. r includes the half-line (Ai, cxj) x {0}. 



Definition 3.2. A subset of the parameter space T C 
if there is some c < oo such that for all (A, a) G F 



^ (2^c/) 



E 

k=2 



dk 



k\ 



-\ak\ < c\X\ 



V ^ is called sub-conical, 



(3.1) 



In this context it is worth noting that under Assumption A2 the interaction obeys 
the bound; 



||W(q:)|| := sup |ZY(x;q;)| < n 



Lk=2 



dk 



k\ 



\ak\ 



(3.2) 



A useful criterion of localization is expressed in terms of the fractional moments of 
the Green function, with the average being carried out over both the disorder and the 
energy within an interval / C M. For this purpose we denote: 



E 



/ ■ 



I/I 



E [.] dE . 



(3.3) 



Definition 3.3. A robust subset of the parameter space, F C M+ x W^^ is said to 
be a domain of uniform n-particle localization if for some s G (0,1) there exists 
^ = £,{s,n,p) < oo and A — A{s,n,p) < oo such that for all (A, a) G F, all 
k G {1, and all x,y G C^"): 



sup sup E/ 

ICR ncz-i 
|/|>i 



Gl^\x,y)\ ) < ^e-'i'"'" (3.4) 
where the energy variable on which G depend was averaged over the intervals /. 
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In the above definition we have incorporated the specific choice of the distance, 
dist-H, only for convenience. As was explained in the introduction, it seems natural to 
expect exponential decay also in terms of the symmetrized distance, distg, but that is 
not proven here. 

Let us also add that localization in the sense of ( |3.4| i implies various other, more in- 
tuitive and physically relevant, expressions of the phenomenon; in particular dynamical 
(Section[4]i as well as spectral localization (Appendix [B). 

In this work we shall focus on the proof of existence of robust regimes of local- 
ization for any finite n and p, without monitoring closely the values of the localization 
length f , and amplitude A. In particular, the subsequent proof yields a localization 
length which degrades heavily when the number of particles n increases. 

Concerning the value of s in the above definition, it is helpful to notice 



Lemma 3.1. If \X\ is bounded away from zero and the condition p.4| i is satisfied for 
some s e (0, 1) then it holds for all other s G (0, 1) at adjusted values of ^ < oo and 
A < oo. 



Proof. Jensen's and Holder's inequality imply that for all r < s < i < 1 



E,[|Gj^(x,y)|nr <E/[|Go(x,y)n 



|Go(x,y)|' 



E,[|Gn(x,y)r 



(3.5) 



The first term in the last line is bounded, E/ 
to (|Z2ll. 



Go(x,y)|*J < Cit)\X\-\ thanks 

□ 



4 Multiparticle eigenfunction correlators and the Green 
function 

4.1 Eigenfunction correlators 

A convenient expression of localization, and also a convenient tool for the analysis, 
is provided by the eigenfunction correlators. By this term we refer to the family of 
kernels, for x, y S (Z'')": 

g[,")(x,y;/;.) (4"' , A^^}(^o"^) ^^^V"^ | (^^ , Ai^}(^o"^) 4"^> f ^ 

£;ecr(Hj'"*)n/ 

(4.1) 

where / C M is a subset of the energy range, C Z"^ is a finite subset, P{E}iH^^) is 
the spectral projection on the eigenspace corresponding to the eigenvalue E, and s G 
[0, 1] is an interpolation parameter This definition extends naturally to also unbounded 
O C Z'', provided ^ffj"' P™^ point spectrum within /. The notation used 

here differs from that of 0Ai94l by allowing for degeneracies in the spectrum. As 
mentioned in the introduction, while the spectrum of a one-particle Hamiltonian with 
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random potential is almost surely non-degenerate, degeneracies do occur in the non- 
interacting multiparticle case. 

When the domain and the value of n are clear from the context, or of no partic- 
ularly importance, the sub/super-scripts on Q may be suppressed. When s is omitted, 
it is understood to take the value s — 1, which for many purposes is the most relevant 
one. 

An essential property of the kernel is the bound (at s = 1): 



sup \{5^, f{Hn)5y)\ < |Qo(x,y;/)| . 

ll/l|oo<l 



(4.2) 



In its dependence on the parameter s the kernel is log-convex, i.e., for any A e [0, 1] 

g(x, y; /; (1 - X)po + Xpi) < Q(x, y; /; so)'^'-^^ Q(x, y; /; si)^ . (4.3) 
Moreover: 

g(x,y;/;0) = 



Q(x,y;/;1) < 



^ (<Sx,P{B}(Fo)<5x) < 1 

E 



\{d^,P{E}{Hn)Sy)\ < 1. 



(4.4) 



where the latter is by the Schwarz inequality. A useful implication of equations ( |4.4[ l 
and ( |4.3| l is that for any < s < i < 1 

Q(x,y;/;t) < Q(x,y;/;s)i^ . (4.5) 

The relations ( |4.4[ ) and ( |4.3| l played a role in the strategy which was used in IIAi94l 
for the deduction of dynamical localization through Green function fractional moment 
bounds. As we shall see next, the method can be extended to many particle systems. 

Most of our analysis will be done in finite volumes. A minor subtlety concerning 
the passage to the infinite volume limit, is that we do not have an a-priori statement 
of convergence in this limit of the eigenfunctions, nor of the eigenfunction correlators. 
Nevertheless, one has the following statement. 

Theorem 4.1. Suppose that the following bound holds for a sequence of finite domains 
f2 which converge to T/^, and a fixed interval / C M, 



E 



Q[r^(x,y;/) 



(4.6) 



with K{-, ■) some kernel, i.e., a two point function defined over the space of pairs of n- 
particle configurations, and some A = A{n) < oo. Then, within the n-particle sector, 
the infinite volume operator H(uj) satisfies: 



E 



sup 

l/IU<i 



(4"), /(i/("))<5(")) 



(4.7) 



Furthermore, if ( |4.6| l holds with K{x.,y) = 2dist-H(x, y)/^ then one may also con- 
clude that the n-particle spectral projection on I is almost surely given by a sum over 
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a collection of rank-one projections on eigenstates which decay exponentially, each 
satisfying a bound of the form: 

|V^(x;^)|2 < A{Lo-n) (1 + |x^|)'"''+' g- dist„(x,x^)/c ^ ^4 8^ 

where A{uj; n) is an amplitude of finite mean, and the decay is from a configuration 
at which the wave function is non-negligible in the sense that 

\^{^,:^)?> ^ (4.9) 



Eye(z^)-(l + |y|) 



With the natural modification, the last statement is valid also in case K{yL,y) is 
given in terms of any of the other distances which were mentioned in the introduction, 
i.e.,dist(x,y) or dists(x, y). 

Except for a minor reformulation of a known bound, this relation is in essence well 
familiar from the theory of one particle localization (it was used already in ||Ai94| ). We 
therefore relegate its proof to the Appendix ( [B] i. 

As it turns out, averages over the disorder of the eigenfuction correlator are closely 
related with Green function's fractional moments. The rest of this section is devoted to 
the relations between the two quantities. 

4.2 Lower bound in terms of Green function's fractional moments 

The following (deterministic) estimate allows to bound fractional moments of Green 
functions in terms of eigenfunction correlators. 

Theorem 4.2. Let C Z''. For any s e (0,1) and any interval / C M 

y^|G[r^(x,y;£;)[ di? < q|,"^ (x, y, M)^ (4.10) 

One may note that this bound is useful only in case of complete localization of 
all eigenfunctions, but that suffices for our purpose. The bound may be improved 
with a restriction of the eigenfunction correlator to a finite, slightly enlarged, interval 
/' D /; the contribution to the Green function from eigenfunctions outside / being 
handled with the help of quasi-analytic cutoff in the sense of Helffer-Sjostrand, and the 
Combes-Thomas Green function estimate. 



Proof of Theorem \4. 2\ We split the contribution to the Green function into two terms 
depending on whether {5^ , P{E}{Hn) Sy) > or ((5x , P{E}{Hn) Sy) < 0, 

G=^(x,y;2):= (4.11) 

Eea(H) 
sign(5^,P|£;,(Hn)5y)=± 
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(Note that the eigenfunctions of may be taken to be real. In the complex case one 
would have four terms instead.) Using \a + fej'* < |a|'* + \b\^ we thus get 

/ \G{^,y;E)\' dE<Y, I \G*{^,t,E)\' dE 

= ^,s/ |{£;e/| |G#(x,y;£;)|>t}| (4.12) 

#=± 

Boole's remarkable formula, which state s that | {x G K | | '^nPni^n ~ > i}| = 
2 Y.n Pn fo"" Xn eM.,Pn,t > 0, IIBo57l . implies that 

|G#(x,y;z)| >t}| = ^ ^ (cJ, , P{s}(i?o) <^y) ^Q#(x, y, M) . 

sign(5,,P{E}(Hn)5y>=# 

(4.13) 

Substituting in integral ( |4.12| i the maximum of ( |4. 13| l and the length of the interval, |/|, 
one arrives at 

J^\G{^,y;E)\' dE < [Q+(x,y,M)^ + Q-(x,y,M)^] 

< Q(x,y,M)^ (4.14) 

1 — s 

□ 

4.3 Upper bound in terms of Green function's fractional moments 

For the proof of our main result we need also a converse bound to ( |4.10| i. 

In the one-particle case there is a simple passage from exponential decay of Green 
function fractional moments to similar bounds on the mean value of the eigenfunction 
correlators, and thus to dynamical localization M Ai94l . In effect, it is based on the 
following relation, which is a somewhat more explicit statement than what is found 
in ||Ai94l . 

Lemma 4.3. For any finite domain C Z'^, x G A, s €E (0, 1) and Borel set / C M, 

f^\x,V-I,.s)\^^^^^^^^^ = \\r' j\^{5,,{H^^'> ~E)-Hy)^ dE (4.15) 



where the left side involves the eigenfunction correlator for the one parameter family 
of operators H^\v) := H^^^ + Xv Px acting in £'^(A). 

The combination of ( |4. 15| l and the 'double sampling bound' ( |2.4[ ) produces the 
desired upper bound on the the expectation of the eigenfunction correlator - in the 
one-particle case. 

To extend this argument to the multiparticle case, an extension is needed of the 



averaging principle which is expressed in Lemma 4.3 Following is a suitable general- 
ization. 
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Lemma 4.4. Let s G [0, 1) and ft C be a finite set, and u a point in fi. Then for all 
X e C(") {VL; u) and y S C^") (^2); 



iV,(x) 



Ke(T(A'„(£;)) 

X |(,5x,n,(£;)v/^(i/j,-£;)-i(5y) 



dE . (4.16) 



where Ii,^{E) = P{k}(-^«(-E')) " f/^^ spectral projection on the eigenspace at eigen- 
value nfor the E- dependent operator: 



K^{E) := ^^{Hn-Ey 



(4.17) 



/(») 



which we take as acting within the range of in Ti^j . 

Since the proof takes one on a technical detour, we have placed it here in the Ap- 
pendix [C] Using this averaging principle, we get: 

Theorem 4.5. Let ^ a be a finite subset, and u € fi. Suppose x,y € C'-"-'(ri) is 
a pair of configurations such that the number -/V„(x) of particles of a at u is at least 
one. Then for any s € (0, 1) and any interval / C M 



iY,(x)Eb[,")(x,y;/,s) 



< 



K\E„ 



\l-s 



E 

ifeC(")(f2:«) 



iV„(x) 



s/2 



E 



dE. (4.18) 



Proof. It follows from ( |2.4| l that the conditional expectation of any non-negative func- 
tion / of the random variables {V{x)}x^n, conditioned on the values of V at sites 
other than u, satisfies: 



K[f{Viu))\{Vix)},^.u] <K\Eo\'E 



(4.19) 



We apply this relation to / the eigenfunction correlator. The quantity which one 
then finds on the right side of ( |4.19[ ) can be rewritten with the help of Lemma |4.4| 
The claimed bound then easily follows using the fact that \a + < |a|* + (for 
< s < l),andj:jS^,U,{E)6^)^-' \(S^,U,{E)S^)\' < hin^JE^. □ 

Applications of the above result are restricted to bounded intervals /. In order to 
control the eigenfunction correlator associated with the tails of the spectrum we also 
use: 

Lemma 4.6. Let ncZ'^andE> 0. Then for every x e C^"^^)- 



E 



< E 



exp{mm{l,{n\X\)-^}{2dn+\\U{a)\\^- E)} . (4.20) 
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Proof. The Chebyshev-type inequality, '^^\[-e,e){x) < e (e + e ), reduces 
the bound to one on the semigroup for which we employ the Feynman-Kac representa- 
tion (cf. IICL90 Prop. 11.3.12]) to show that for any t > 0: 



(4.21) 



< 



J E exp ^ (^Y. AF(u)iV„(y(s)) +Z^(y(s);a)j ds 



,(x;t) 



(dy) 



< E 



,nt\\\\V{0)\ 



exp{t{2dn+\\Uia)\\J} (s^ , e'^^ "^^'-^ 6, 



where ^'q''*' is the measure generated by J2j ^fi.j on path {y(s)}o<s<t> starting in and 
returning to x in time t. The last inequality is a version of Jensen applied to average 
(nt)^^ SuGJi /o (■) ^!i(y(*)) ds in the exponential. A similar bound holds for t < 0. 
The proof is completed using the operator bound (^x , e*'^"(5x) < e^''"'*' for the free 
semigroup and the choice t = min{l. (n|A|)~^}. □ 



5 Implications of localization in subsystems 

In the induction step of the proof of the main result, we shall be considering for a 
system of n particles the consequences of localization bounds which are akeady es- 
tablished for subsystems. In this section we present some results which will be useful 
for that purpose; first considering the case when the subsystems are combined without 
interaction, and then the more involved situation where the two subsystems are cou- 
pled via short range interaction. For a partition of the index set {1, n} into disjoint 
subsets J and K, we denote the coordinates of the two subsystems as xj := {xj}j^j 
and correspondingly xk- 



5.1 Localization for non-interacting systems 

When two subsystems are put together with no interaction, the Hamiltonian is - in 
natural notation, 

<-^=<©i//f\ (5.1) 

acting in £^(ri)l-^l ^^(fi)!^'. A complete set of eigenfunctions of the operator sum 
can be obtained by taking products of the subsystems' eigenfunctions. Clearly, if the 
subsystems exhibited spectral locahzation, that property will be inherited by the com- 
posite system. 

The question of localization properties of the corresponding Green function, which 
is an important tool for our analysis, is a bit less immediate: g[/'^' is a convo- 
lution, with respect to the energy, of the subsystems' Green functions, i.e., for any 
X = (xj, xk), y = (yj, yx), and z e C\R: 

G[/^^^(x,y;z)= / G[^\^j,yj;z - OG^nH^K,yK;0 , (5.2) 
Jc 
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where C is any closed contour in C, which encloses the spectrum of but none 

of H^^^ — z. Given the singular nature of the i?-dependence of the Green function, 
localization in the sense of Definition 3.3 is not immediately obvious. To establish it, 
we take a detour via eigenfunction correlators. These are less singular in E, but share 
the convolution structure, which in this case can be written in the form: 



Q[/'^)(x,y;/)^ !( 



< q[/^ (xj , y J ; M) Q[f ) (x^f , y K ; M) , 
The following result will allow us to apply this relation. 



Lemma 5.1. Let T e 



(5.3) 



X M.P^^ be a sub-conical domain of uniform n-particle 
localization. Then there exist A, ^ G (0, oo) such that the eigenfunction correlator 
corresponding to up to n particles, i.e., k g {1, . . . , n}, is exponentially bounded for 
all (A, a) G r, and all x, y G {Z'^)'' : 



sup E 



(5.4) 



Proof. As an immediate consequence of ( |43| l. Theorem 4.5 and Lemma A. 3 we know 
that there is A, ^ G (0, oo) such that 



E 



0[f)(x,y;[-i?,£;]) 



distH(x,y)/5 



(5.5) 



for any E > Q and any (A, a) G F. For a bound which is uniform in E we combine 
this with Lemma 4.6 which with the help of the Cauchy-Schwarz inequality implies: 

E 



< [E 

< E 



E 



Q[';\^,y;[~E,E]) 



Q[';\^,^;R\[-E,E])\ E \Q',';\y,y;R\[-E,E]) 



1/2 



exp {min{l, {k\\\)-'^} {2dk + \\U{a)\\^ - E)} 



(5.6) 



Choosing the cutoff at _E = 2dk + \\hl{a)\\^ + max{l, A: |A|} dist-H(x, y)/^ one 
obtains the claimed exponential bound. Using the fact that F is sub-conical in sense of 



Definition 3.2 the above argument yields a (A, a) -independent amplitude "yl" in ( |5.4| l. 

□ 

The two-way relation between the eigenfunction correlators and fractional mo- 
ments of the Green function, and the factorization property ( |5.3| l, allows us now to 
establish: 

Theorem 5.2. Let F G IR+ x Mp~^ be a sub-conical domain of uniform n-particle 
localization. For any s G (0, 1) there is ^, A G (0, oo) such that the Green function 
of the composition \5.l) of any pair {J,K) of systems of at most n particles (i.e.. 
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max{| J|, jii'l} < n) is bounded for all (A, a) G T and all x = (xj,xx), y = 



sup sup E,f G[/'^)(x,y) < Ae-''-^<'"^--y^/^ , 
|/|>i 

M'!f/idisty'^^(x,y) niax{dist-H(x,/,yj),dist-H(xif,y/f)}. 



(5.7) 



Proof. By Theorem 4.2 and the Jensen's inequality, for any s £ (0, 1) there is a 

C 



constant C = C{s) < oo such that 



G 



(J,K)) 

n 



(x,y) 



(x,y;/) 



(5.8) 



The claim follows by combining: /) the product formula ( |5.3| l, //) the uniform bound 
Q < I, and the bound of Lemma 5.1 applied to the factor with the greater separa- 
tion. □ 



5.2 Decay away from clustered configurations 

We now turn to the more involved situation, where a system consists of subsystems 
which separately exhibit localization, but which are put together with an interaction. 
Intuitively, the decay of the fractional moments for the subsystems should imply small- 
ness of the corresponding kernel for the composite system for pairs of configurations 
where at least one of the pair can be split into two well separated parts. 

To express this idea in a bound, we shall use the notion of the splitting width of a 
configuration: 

^(x) max min Ix, — Xfel, (5.9) 

J.K jeJ,keK ■' 

JOK={l,...,n} 

where the maximum runs over all the two-set partitions of the index set {1, . . . , n}. It 
is easy see that diam(x)/(n — 1) < ^(x) < diam(x) (cf. Appendix [a|. 

Theorem 5.3. Let T^^}_^ C M+ x M.p^^ be a sub-conical domain of uniform {n — 1)- 
particle localization, with n > 2. Then there exist some s € (0, 1), A, ^ < oo such that 
for all (A, a) € r|f2i andallyi,y e {I'^Y : 



sup sup E/ 

/cR ndz-i 
|/|>i - 



|G[r^(x,y)r 



< Aexp(--min{distT^(x,y),max{^(x),£(y)}}) . (5.10) 



Proof We fix X, y G C(")(r2) and assume without loss of generality that ^(x) > ^(y). 
We then split x into two clusters, xj, xj^ such that 



£(x) 



mm a;,- — X]^\ 
jeJ.keK 



(5.11) 
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Between the two clusters, xj, xk we remove all interactions, leaving the inter-cluster 
interaction untouched. The resulting operator is a direct sum of two non-interacting 
subsystems, H^^'^^ := iJ^^' ® H^f^ acting in fi^^'^^ ® where 

\J\ \J\ 
i^^-'^^JI-^J E Ua{MaM) (5.12) 

i = l fe=2 ACZ'*:|A|=fe 

diam A<£[/ 

and similarly for H^^\ Denoting the Green function corresponding to Hq'^^ by 
Gq '^\ and using \a + b\'^ < |a|'* + and we thus have 



E, 



G[/'^)(x,y)-Gn(x,y) 



(5.13) 



Since G^^^'^^ is a Green function of a composite system, whose parts are assumed to 
exhibit uniform (?i — 1) -particle localization. Lemma 5.2 guarantees the existence of 
s e (0, 1) and A,£, e (0, oo) such that for all (A, a.) eT^\ and all x, y: 



Er 



Gr^(x,y) 



< Ae' 



■ dist. 



^'■'^'(x,y)/? < ^g-dist„(x,y)/? 



where the last step is by the general relation dist^'^'' (x, y ) > dist-^ (x, y) . To bound 
the second term we use the resolvent identity 

A:= G[2''^)(x,y;z)-Go(x,y;z)= ^ G^^'-''\x,w; z) Uj^^) Gn{w,y; z) , 

wGC(")(a) 

(5.15) 

where L/j.A' — H —H'-'^'^\ In order to be able to apply the Cauchy-Schwarz inequahty 
we first decrease the exponent with the help of Holder's inequality, 

l+2a 1-s 

(e/[|A|^ 



(5.14) 



E/[|A|''] < (E/[|A|5fTTifT 



]) (E,[|A|5]j <c\X\-^i^ 



(5.16) 

where the last inequality is due to ( |2.2| i. Inserting ( |5.15| l and using the Cauchy-Schwarz 
inequality together with (|2.2[) we thus obtain 



(5.17) 



E,[|Ar] < ^ J2 Pj.Ki^r^ (e [ G«^)(x,w;z) 
where we abbreviated /3 := 2(2+ s) • ^'^ estimate the right side note that 



sup|f/j,if(x)| =sup|^afc ^ C/a(A/a(x)) 

fe=2 ACl.'^:\A\=k 
diam A<£u 

X 1 [There is j £ J,k £ K s.t. \xj — Xk\ < £u] 



< \\U{a) 



(5.18) 
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Moreover, the distance of x to the support of Uj^k is bounded from below, 



inf dist^^^-^' (x, w) > €(x) - £u : 

wSsupp U.I^K 



(5.19) 



by the triangle inequality and ( |5.11| l. As a consequence, ( |5.14| i and Lemma [A3] yield 



E/[iAn 



|A| 



\\ll{a)\ 



exp <{ ^^(x) 



< 



A exp|-^disty^^^(x,w)| <(X3. (5.20) 



The proof is concluded by noting that ||Zi(Q:)||^ < C|A| for some C < oo for all 
(A, a) g r|f2i since the latter is sub-conical. □ 

The above result will allow us to insert in sums over n-particle configurations a re- 
striction to ones of a limited diameter. The following bound will be useful for estimates 
of the remainder. 



Corollary 5.4. Under the hypothesis of Theorem 5.3 there exist s G (0, 1) and A,£^< 
oo such that for all < r' < r, the quantity 



ksi^,r,r') := sup 

\x — y\>2r 



^ E[|Go(x,y;z)n (5.21) 



yeC<"'(n:y) 
eC<"'(0;a;)\C<';'(n;a;) 



satisfies for all (A, o;) e ^lili- 

ks(Sl,r,r') < Ar-''("-i' |f^|""iexp 



(n-l)e 



Proof By Lemma A. 2 for all x, y in the sum in ( |5.21| i: dist-H (x, y) > |x — yj — r > 
r > r'/(?i — 1). Theorem 5.3 hence guarantees that for some s e (0,1) and A,£_ < oo: 



E[|Go(x,y;z)r] < Aexp 



(5.22) 



The proof is completed by bounding the number of configurations in ci^\fl; y) and 
C(")(f^;2;) by \Crln;y)\ < n(4r)''("-i) and \C{n;x) \ < n respectively. □ 



6 Proof of the main result 

We now turn to the proof of Theorem |1.1| That is, we shall prove that there exists 
a monotone sequence of decreasing sub-conical non-empty domains in the parameter 
space, 

ri^^ 3 ••• 3r,[,f) D ... , (6.1) 
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such that for each n, the Hamiltonian with parameters in tI^^ exhibits uniform 
n-particle locaHzation in the sense of Definition |3.3t each of the domains including 
regimes of large disorder and of weak interaction. The proof will proceed by induction 

on n, the induction step consisting of a constructive restriction of the domain. 

(v) 

Establishing a sub-conical domain r„ of uniform n-particle localization would be 



sufficient for our purpose, since Lemma 5.1 then implies that there is A, ^ e (0, oo) 



such that for all (A, a) e vi^K all A: e {1, . . . , n}, and all x, y G {Z'^)'': 

sup Q|^''''(x,y;M) < A e" '^''''"(^'^^/^ . (6.2) 



Spectral and dynamical localization for H^''\ as claimed in 1. and 2. of Theorem |l.l[ 



follow using Theorem 4.1 Assertion 3. on the shape of the nested decreasing domains 



of uniform localization will be verified in the course of the inductive construction. 



6.1 Analyzing clustered configurations 

An essential component of the proof is to show the exponential decay of the finite 
volume quantities: 



Bi^^\L) sup sup |9Al| V V E, G^^\^,y) 
|/|>i " ye^At xec^2 (f^'O) 

yeC(">(n;!y) 



(6.3) 



where A^ := [— L, LY H Z'*, and s e (0, 1). One may note that the sum is resticted 
to configurations in the form of separate clouds of particles with diameter less than 
:— L/2, which are guaranteed to include at least a pair of sites at distance L G M 
apart. By the Wegner estimate ( |2.2| l, for all s G (0,1) and L e N: 



where the constant C = C{s, d) < oo is independet of (A, a) £ M". 

The following rescaling principle will be used to show that B'f'\L) decays ex- 
ponentially provided that there is some Lq for which it is sufficiently small. In its 
formulation, we consider length scales which grow as 

Lk+i := 2(Lfc + 1) , (6.5) 

i.e., Lk = 2'' {La + 2) - 2 for all k e Nq. 

Theorem 6.1. Let t[^}_^ C M-|_ x M^^^ be a sub-conical domain of uniform [n — 1)- 
particle localization, with n>2. Then there exists s G (0, 1), a, A,p < oo, and > 
such that 

Bi-HLk+i) < ^ Bi-^L^r + ALf^, e-''^^^ (6.6) 
for all (A, a) e r|f2i and all k G Nq. 



21 



In order to keep the flow of the main argument clear, we postpone the proof of this 
assertion to Subsection l6.3l 

Quantities satisfying rescaling inequalities as in ( |6.6| l are exponentially decreasing 
provided that they are small on some scale. This is the content of the following lemma. 
For its application we note that 

S{Lk) := Bi'^HLk) , Lk := 2\La + 2) , (6.7) 

satisfy ( |6.8| l. 

Lemma 6.2. Let S{L) be a non-negative sequence satisfying: 

S{2L) < a S{Lf + h L^P e^^"^ (6.8) 

for some a,b,p [0, oo) and ly > 0. If for some Lq > there exists 77 < oo such that 
2P 

1. rf' > ab + rj 



TP ' 

^0 



2. 1 > aS'(Lo) +?7Lge 



then for all k G Nq.' 

S{2''Lo) < fl-i exp (-/i2'=Lo) . (6.9) 
Proof From it follows that the quantity R{L) := aS{L) + r]LPe~''^ satisfies 

R{2L) < {aS{L)f + (<^b + rj'^^ L^Pe"^"'' 

< {aS{L)f + Tj^L^Pe-^-"^ < R{Lf . (6.10) 
The claimed ( |6.8| l follows by iteration, using 

6.2 The inductive proof 



Proof of Theorem 1.1 As explained above in this section, it suffices to establish sub- 
conical domains of uniform n-particle localization (in the sense of Definition |3.3[ l. 

As an induction anchor we use the fact 0AM93II that there is some Ai e (0,cx)) 
such that r^''^ := (Ai , 00) X ]RP ^ serves as a domain of uniform localization for the 

one-particle Hamiltonian ij'^^ (Note that the last {p — l)-components of T^f^ are 
irrelevant for H^^\) 

In the induction step (n — 1 — > n), we will first construct a robust, sub-conical 

(v) (v) 

domain r„ ^ C r„_j^ and pick G N sufficiently large such that at some inverse 
localization length /i > 0: 



for all k G No and all (A, ot) e tI{'\ Based on the induction hypothesis. Theorem 
and Lemma [672] this is will done separately in the two regimes of interest: 



6.1 
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Strong disorder regime: Here we choose G N by the condition 2^6 < 1/4. 
Based on that we pick rj in the range 

2P+1 1 r 

-Ti^<V<WYpe (6-12) 

which is non-empty by our previous choice of Lq. Notice that for this choice 
r]2P/L^ < 77^/2. We now restrict the domain T^li by choosing A„ > A„_i 
large enough such that for all A > A„: 

and ^Bi-HLoXl (6.13) 



This is possible by ( |6.4| l and the fact that A is independent of A. Lemma 6.2 
hence guarantees that 

log(aA--i3i"^(Lo) + l/2) ^ log(2aA--i3^)(£o)) ^ ^ ^^^^^ 

Lq 2Lo 

serves as an inverse localization length in ( |6.11| l, which is valid for all (A, ct) in 
the following regime of strong disorder: 

Dir'' := {(A, a) I A > A„} n T^fli • (6.15) 

Weak interact ion r egime: The localization bounds in MAM93II and arguments as in 



Theorem 



5.2 



ensure that _Bi"''(L) ^ as L — > 00 for a = and all A > Ai. 
We then pick Lq E N and ij such that for a = 0: 

^B(")(Lo) + ??ige-''^°<l, (6.16) 

-^+Vjp<V^ (6.17) 

Since the finite-volume quantity B^s^\Lq) is continuous in ol at 0, for every 
A > Ai there exist q;(A)j > 0, j € {1, . . . ,p} such that and ( |6.17| l is 

maintained for all (A, a) in the following regime of weak interaction 

:= {(A, a) I A > Ai and [a^ l < a(A)j for all j € {1, . . . ,p}} n r|f2i • 

(6.18) 

By Lemma 6.2 /j, as in ( |6.14| l with A„ replaced by Ai hence serves as an inverse 
locahzation length in ( |6.1 l| l, which holds for all (A, a) £ 

Summarizing, we have thus established ( |6.1 l| l for all (A, a) e I?!^^ U /^^^ and a 
suitably large Lq E N. By further restriction, we may find thus find a robust, but 
sub-conical domain 

(0 ^) r(f ) C D(f ) U /t^') (6.19) 

for which ( |6.1 l| l holds. 
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To complete the induction step, it requires to establish the exponential decay p.4[ l 
of the n-particle Green function for all (A, a) E 

For this purpose, we select x E X and y E Y such that dist->^(x,y) = |x — ?/|. 
In case |x — > Lo, which we may assume without loss of generality (wlog), there 
exists a unique fc G No such that the box A^^ (x) which is centered at x satisfies: 

• y ^ ^Lk{x) and y E Al 

fc+i l-^/' 

• Lk < \x — y\ < cLk for some c E (0, oo). 

We may furthermore assume wlog that diamx < Lk/2 is sufficiently small such that 



X E C'"-'(ri n Ai^,(a;)), since otherwise p.4| i (with k = n) follows from Theorem 5.3 
The resolvent identity, in which we remove all the terms in the Laplacian which connect 
n A^^ (x) and its complement, then implies 



E,[|Ga(x,y)r]< E 

wec<"'(r!nAt ;«) 



X Er 



c 



^ ml ^ _^^P E E 



[I I ^ 

|GonAi^^(x)(x,w)| |Go(w',y)|'' 

E, [|Gji(x,w) 



(6.20) 



where i?w,w' are the matrix elements of the Hamiltonian between and 6^'. In- 
equality ( |6.20| i follows from p.2| l by first conditioning on all random variables apart 
from those associated with sites y and u' which are both outside A^^. (x). The sum in 
( |6.20| l is split into two parts, depending on the diameter of the configuration w: 



C ^ 
E 



E 



Er 



4"^(x,w; 



(6.21) 



By ( |6.11| l both, in the strong disorder regime and the regime of weak interactions, the 
first term on the right side is exponentially decaying in L^. To bound the second term 



we employ Theorem 53] again, 
E/ 



G[r^(x,w) 



< A exp 



2(n-lK 



(6.22) 



since dist-H(x, w) > \x ~ u\ ~ diam(x) > Lk/2 and ^(w) > Lh/{2{n — 1)). The 
sums in the second term have at most |9Alj.(x)| and n|Ai^(a;)|""-'^ terms, respec- 
tively. Hence, this second term is also exponentially decaying in Lk and hence in the 
Hausdorff distance, dist^^ (x, y). This completes the proof of the induction step. □ 
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6.3 Proof of the rescaling inequality 



In the above proof we have postponed the derivation of Theorem 6. 1 which provides an 
essential step for estabHshing regimes of exponential decay. We shall do that now. The 
construction and the analysis are inspired by arguments familiar from the one-particle 
setup. The following lemma essentially extends arguments in IIAS+Oll . 

Lemma 6.3. Let n C Z'^ and V,W C n with dist(y, W) > 2. Then for all x € 



E 



4"^(x,y;^) 



< 



C 



Gv^(x,w;z) 



Gl^^(v,y;z) 



(6.23) 



■ E E E 

uGdV\an vreC'-"^ {V;u) 

X E E ^ 

vGdw\anvec("HW;v) 

where the constant C — C{s,d) < oo is independent of (\, a) G MP. 

Proof. Let H^^i denote the matrix element of iJ^ between 5^ and . A twofold 
application of the resolvent identity, in which we remove all terms in the Laplacian 
entering which connect sites in V with ^l\V and W with ^l\W , yields 

Go(x,y;z)= ^ Gy(x, w; z) iJw,w' Ga(w', v'; z) i/v,v' Gvi/(v, y; 2) . 

wec<"'(v),w'ec<"'(n;r!\y) 
vgc<"' (w),v'ec<"' (0;0\w) 

(6.24) 

Using |a+6|^ < |a|'' + jfol" for any s G (0, 1), we proceed by establishing the following 
implication of Theorem |2.1[ 

E[|Gy(x,w;z)r |Go(w',v';z)r |G^(v,y;z)r] 

< ^E[|Gy(x,w;z)n E[|GH'(v,y;z)n . (6.25) 

For its proof we note that there are two sites u' ,v' ^ V U W for which the configu- 
rations w' and v' in the sum in ( |6.24| i have a particle at u' respectively v'. We may 
therefore first condition on all random variables {V{x)}^^^^, and use ( |2.2| i to esti- 
mate the factor in the middle. The other factors are independent of each other and of 
V{u'), V{v'), such that the expectation value factorizes. 

The proof of ( |6.23| l is completed by noting that for any fixed w the number of 
w' G C(") {n; n\V) with Hw.w' 7^ is at most 2d. □ 



We are now ready to give a 



Proof of Theorem \6. 1 \ We first restrict the sum in the definition of B^J^\Lk+i) to con- 
figurations X, y which have a diameter less than r^^. 



laAi.+J sup 

|/|>1 



sup 



E E 



Er 



-(") 



G[r'(x,y) 



(6.26) 
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The error is controlled with the help of Corollary 5.4 since y E dKLk+i we have 
> Lk+i and hence 



+1 

2(in-2 
fc+1 



exp 



2(n-lK 



(6.27) 



The sum B'f'\Lk+i) is now estimated with the help of Lemma 6.3 in which we pick 
y = Alj. n and W — A^^ {y) n 51, where the last box is centered at y G dhLk^ri ■ 
Since (fi; 0) C cif^ {V; 0) and cif (fl; y) C C^") (M^; y) we thus obtain 



Bl"^(ifc+i) < 



9 A 



sup 

/CR 

I^I>1 



sup 



E E 



Er 



I - uedAL^ xec (n;0) 

weC<"'(f2;u) 



G[r^(x,w) 



X sup sup E E 

vGCf"'(n;D) 



4"^(v,y) 



(6.28) 



Thanks to translation invariance we may shift y to the origin in the last line. Moreover, 
|(9Ai^.^j I < 4^*^^ I^Alj. I and we may again restrict the summation to configurations 
with a smaller diameter using Corollary |5 .4| again. 



2d-2 



< 2 



4d-3 



c 

W 
c 



S(")(Lfe) + |aALj2fc,(Ai,,ri,,riJ 

Lk 



(6.29) 

□ 



Appendix 

A Some distances and separation lemmata 

Following are some natural lengths associated with n-particle configurations, and some 
elementary geometric estimates which are of use within this work. 

In general, for a configuration which is denoted by a bold lower-case letter, we shall 
use the corresponding capital letter to denote its footprint in Z'^, which is the subset 
which it covers; e.g., for x = {xi, ...,a;„} C (Z**)" we let X = UiLil^^i} ^ Z''. 
Also, for subsets the index set, ,/ C {1, .., n}, we let Xj :— IJigjI^^i} ^ (Z'^). 
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A.l Splitting width 

Two convenient measures of the spread of a configuration x = {a;i, a;„}, are: 

1. the diameter, diam(x) := maxj^ke{i,...,n} |xj — xk| 

2. the (maximal) splitting width, which we define as the supremum over r for which 
there exists a partition of the set X into two subsets at distance r apart, or: 

ehc) := max dist(Xj, Xr) , (A.l) 

J,K: JUK={l,..,n} 

where dist(A, B) = minygA, veB |u — v|, with \u — v\ the Euclidean distance. 

Since both quantities depend only on the footprint of x, in a harmless abuse of notation 
we may also refer to diam(x) as diam(X) and to £(x) as 

Lemma A.l. For any configuration x e (Z'^)".- 

^^-^diam(x) < ^(x) < diam(x) . (A.2) 

Proof. The upper bound is totally elementary. To prove the lower bound on £(x) con- 
sider the one-parameter family of sets 

n 

Xr ■■= [jiy&I^" : \y-xj\<r}. (A.3) 

For any r > such that is connected one clearly has 

max \xj — Xk\ < 2r(n — 1) . (A.4) 
j,fce{i,...,n}' 

It follows that for any r such that 2r < diam(x) /(n — 1) the set Xr is not connected, 
and hence there is a partition of the configuration x into two subsets whose points are 
at distances greater than 2r. This implies that also 2r < £(x). Optimizing over such r 
we find that diam(x) /(n - 1) < f (x). □ 



A.2 Distances in tlie configuration space 

In addition to the regular distance between subsets of which is mentioned above, 
there exists also the notion of the Hausdorff distance, which is defined as: 

dist7^(X, Y) := max{maxdist({u}, Y), maxdist({v}, X)} , (A.5) 

u£X vGY 

for any X, y e 1'^. 

In a slight abuse of notation we shall employ the symbol dist^^ also for the induced 
Hausdorff distance between configurations (i.e., x, y e (Z'*)"): 

distH(x, y) := distH(X, Y) . (A.6) 
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This distance is clearly most sensitive to the outliers. Another notion is the sym- 
metrized distance: 

n 

dists(x,y) min V |xj - y^j | , (A.7) 

with Sn is the permutation group of the n elements {1, n}. 

The following is an elementary consequence of the definitions. 

Lemma A.2. Let 51 C Z'' and u,v £ ^. For any two configurations x e C*-"-* (fi; u) 
and y g C^"'^(f2; v) , which have a particle at u and, respectively, v: 

dist-H(x,y) > max{dist({u}, F), dist({w}, X)} > |m— wj— min{diam(x), diam(y)} . 

(A.8) 

For convenience let us also place here the bound: 
Lemma A.3. Let C Z'' one/ x e (Z**)". 

7. For any site u £ CI with dist({u}, X) = L, and any ^ > 0." 

J2 e-'^'"'"^^'^)/* < Cmax{L, ^}'*("-i)e--^/«. (A.9) 

yeC(")(0;ti) 

with a constant C — C{n) < oo. 
2. For any ^ > 0." 

g-dist„(x,y)/C < (^^^Q^ 

yec(")(n) 
/or sot7ie C = C{n, d) < oo. 

Proof. It is convenient to use the equality: 

J2 e-^'^t^-^y^/i ^ ^ l[dist«(x,y) < r]. 

yeC(")(n;«) ^ yeC<")(0;u) 

(A. 11) 

The estimate the sum on the right, we note that the configuration y needs to have one 
of its points at u, and the rest (n — 1) points are all within the distance r from X. A 
simple estimate yields: 

i: *M„(x..) <,]<{" III ,A.,2) 

where 2 could also be replaced by bd, which is the maximal value of b such that any 
sphere in M.'^ of radius r includes not more than br'^ lattice points (of Z"^). Substitut- 
ing ( |A.12[ ) in ( |A.ll| i one readily obtains the claimed bound ( |A.9| l. The second claim 
follows from the first by summation over u. □ 
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B From eigenfunction correlators to dynamical and spec- 
tral information 



In Theorem 4. 1 we presented a known method MAi94l for the derivation of information 
on the dynamical and spectral properties of an infinite-volume operator from bounds 
on the eigenfuction correlators of its restrictions to finite domains. For convenience, 
following is an outline of a proof of this result. 

We recall that the assumption is that for a sequence of finite domains which 
converge to Z"^, and a fixed interval /: 



E g[,"^(x,y;/) 



< ^ 



-^f(x,y) 



(B.l) 



with some kernel K{x, y) (e.g., if (x, y) ~ dist(x, y)/0- 



Proof of Theorem 4.1 i) Through a trivial extension of the finite volume operators, 
they can be naturally viewed as acting in the same space £^(Z"'') (acting as on func- 
tions supported outside 51"). Using the Combes-Thomas estimate on the Green func- 
tion, one may see that the operator H{lu) is the limit, in the strong resolvent sense, of 
any sequence of Hn{uj), as ^ Z"* (allowing in the process also arbitrary self-adjoint 
boundary conditions at the receding boundary). Thus, a bound like ( |4.7[ ) but modified 
through a restriction of / to continuous functions can be deduced from eq. ( |4.2| l and the 
general properties of strong resolvent convergence (cf. IIRS79I Thm. VIII. 20]). Since, 
by the Wegner estimate, the mean density of states is a continuous measure, Lusin's 
approximation theorem allows to extend the resulting bound to all measurable and 
bounded functions, thus yielding ( |4.7| i. Of particular interest is the implied dynamical 
localization bound: 



E 



sup \{6^, Pi (H) 



-itH tr 

e Ov 



(B.2) 



ii) By the RAGE criterion (see, e.g. MRS79I ) the projection on the continuous spec- 
trum in the interval / satisfies: 



K[\\Pj..,ont{H)d^ 



= E 



lim lim ^ r J2 K'^- ^^(^) '^y) r '^^ 



yidist (x,y)>_R 



< lim 

R — ^oo 



E 



E 



y:dist(x,y)>fi 



sup \{d^, Pi iH)e 



-itH 



(B.3) 



where the last inequahty is by Fatou's lemma and the natural bound (1) on the summed 
quantity. Thus, under the assumption (|4.6|l. 



E[\\Pi-cont{H)6^\\^] < lim V Ae 

y:dist (x,y)>fl 



A:(x.y) 



(B.4) 
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In case i^(x, y) — 2 dist^ (x, y) / ^, or any of the other distances discussed here (which 
are only larger), the above limit vanishes. Since {(5x}xe(Z'')" is a spanning collection 
of vectors, one may conclude that, under the assumption ( |4.6| l, within the n-particle 
sector H{uj) has almost surely no continuous spectrum in the interval /. 

Hi) Under the above assumption, we will construct a complete set of exponentially 
bounded eigenfunctions which form a subset of {P^e}{H{uj)) 5y\EE (j{H{uj)), y S 
(Z'')"}. Clearly functions in this collection are either zero or eigenfunctions of HioS). 

For the complete set we chose functions corresponding to configurations y € (Z'^)" 
which are E-representative in the sense that 

Exe(z^)"(l + |x|) 

Note that for any E S a{H{uj)) there is at least one i?-representative configuration, 
since J2y{^y ' P{e}{H{uj)) Sy) > 1. We claim that the corresponding collection of 
normalized eigenfunctions: 

V'iJ.y := with E e a{H{u;)) HI , y e (Z"^)" ^-representative, 

(B.6) 

spans the full subspace of eigenfunctions of H{uj) corresponding to eigenvalues in /. 
For if not then there exists E G a{H{uj)) n / and a normalized function satisfying </> = 
P[E}{H{uj))(j>, which is within the orthogonal complement of the subspace spanned 
by ( |B.6| l. This would imply the following contradiction: 

1 = J2 \{^,P{E}{H{u^))Sy)\' < J2 {Sy,P{E}iHi^))Sy) <l. 

ye(z'')"not ye(z'')"not 
iJ-representative i5-representative 

(B.7) 

For bounds on the eigenfunctions in ( |B.6| l, one may apply the Wiener criterion 
which (combined with the Fatou lemma) yields 



E 



J2 \{S^,PiE}{H)Sy)\' 

Eea{H}r\I 

^\{S.,Pi{H) 



< lim E 



< 



' dt 



< Ae-'^'"*«(^'y)/^ (B.8) 



as a consequence of ( |B.2| i in case K{x, y) = 2 dist-H(x, y) /^. Summing the resulting 
bound we get: 



E 



E 



„dist„(x,y)/C 



EeaiH)n 
x.yeCZ")-' 



,(i + |y|) 



nd+1 



\{6^,P{E}{H)6y) 



< oo . 



(B.9) 
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Thus, using the Chebyshev principle, there exists a positive function A{oj; n) of finite 
mean such that for all E G a{H{uj)) n / and x e (Z'')": 

IV'ij.ylx)!^ l[y e (Z'')" is ^-representative] 

< A(tj;n)(l + |y|)^"'^+^e-'^'''*«('''y)/«. (B.IO) 

Since for any E E a{H{uj)) and any y in the above collection the function V'B.y is 
non-negligble at y in the sense of ( |4.9| ) this proves the last claim which is made in 
Theorem O □ 



C An averaging principle 



In the proof of Theorem |4.5| we made use of an averaging principle, which is useful for 
conditional averages where the value of the potential at a single site, u £ Z"*, is redrawn 
at fixed values of the other (random) parameters. This provides a generalization of 
Lemma [43] which is suitable for multiparticle systems. Following is its derivation. 
In the statement of the result, use is made of the holomorphic family of operators 



Kuiz) = VNuiHn-z)- 



(C.l) 



(n) 

which we take as acting in the range of within TY^j . The operator valued function 
is analytic in C\<T{Hn)- For real z = E £ R the operators are self adjoint, and for 
convenience of the (local) argument which follows, we employ an auxiliary index v 
to label the eigenvalues (k^{E)), counted without multiplicity, and the corresponding 
projection operators (XI^^ (_E)). Questions of order do not matter here since we shall 
always be summing over v. 

Along M\cr(i/n) analyticity implies that for all but possibly finitely many excep- 
tional values of E, at which level-crossings occur, both eigenvalues and projections 
may be analytically continued to a small neighborhood of E on which the spectral 
representation 

Ku{z) = ^K,(z)n,„(z), (C.2) 

holds. The operators H^,, (2) are projections, satisfying nK,,(2;) n„_^, (z) — X\K.^{z)5vy , 
though they are orthogonal projections only for real z; cf. IIKa95l Ch. II]. 



Lemma C.l ( = Lemma 4.4 ). Let s e [0, 1) and il C Z'^ be a finite set, and u a point 
in n. Then for all yi£ C'^Q.-u) andy G C(")(r2).- 



dv 

V(x)^V{x)+v \v\^ 



7V„(x)i+4 / Qo(x,y;/,s) 



|A|i 



(C.3) 

dE. 
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Proof. Let us consider the family of operators Hn{v) := + XvNu with the extra 
parameter v, which in effect modifies the value of the potential Vu of Hq- A standard 
resolvent identity leads to the Krein formula: 

y/K[Hniv)-z^ ^ ^{\ + \vKu{z))-^ y/N~u{Hn~z)-\ (C.4) 

For each w, the family of operators \/Nu (^H^{v) — z^ is an analytic in z G 

C\(T(iJn(w)), with residues given by \/Nu P{ b}, where P{e} are the projection oper- 
ators on the eigenspaces of ifo(u), at eigenvalues £' e (T{Ha{v)). The operator valued 
function on the right-hand side of ( |C.4[ ), is singular if and only if either: 

- (Af)-i e <J{K^{E)) (C.5) 



or 



E € a{Hn) 



(C.6) 



and we argue next that the singularities at a{HQ) are removable for almost every value 

of u e M. 

Since the spectrum of Hq{v) is monotone non-decreasing in v, its spectral projec- 



tions can be decomposed (P{ = P{e}'^' ~^ ^{e}^ ^ ^ '^{Hu>{v))) into a part for 
which the corresponding spectrum is strictly monotone and another corresponding to 
spectrum which does not move with v. (Monotonicity plays here only an auxiliary role, 
and could also be replaced by analyticity in v or just smoothness.) At Lebesgue almost 
every w € M the spectrum of Hq{v) corresponding to -Pj^j"'' is disjoint from a{H(i), 
in which case the singularities which {Hi-i{v)—z)^^ may have at a{Hn) are only due to 

the fixed part of the spectrum, with the corresponding residues of -\/iV„ {Hn{v) — z)^^ 
5(fix 



(fix) 



being given by y/ Nu ■ However, functions in the range of the projections P^^^') 



are annihilated by Nu, and thus ^/N^ P^^^' — 0. Since \/Nu {Hn{v) — z)~^ has only 
simple pole singularities, the vanishing of the residue at E implies that the singularity 
of the expression on the right-hand side of ( |C.4| i is removable there. 

In conclusion, for almost every v E M, even if there is an overlap in the spectra of 
Hn{v) and Hn, ^/N:^P{e} = for all E e <T{Hn). 

We now consider E ^ <t{Hq) at which ( |C.5| l holds. Such energies will not coincide 
with exceptional points of level-crossing for Ku{E) for almost all w e M. Therefore a 
simple residue calculation based on \C2\ yields: 



'Nu P{E} = ^VT^ 



Xz)^u{Hn-E)- 



(C.7) 



where k^{E) is that eigenvalue of Ku{E) at which — (Aw)^^ = k^{E) holds, and 
k'^{E) is the derivative of the eigenvalue with respect to E evaluated at that particular 
point. In particular. 



^N,P, 



{E}^ 



'Nu 



d 

Je 



K,{E)-^ 



li^AE). (C.8) 
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Using the relation 5{ g{E) ) g'{E) = J2ueg-'^{{o}) ^(-^ ^ ^^'^ ^^^^ f*-"" ™y 

/ which is continuous on a neighborhood of a{HQ): 

J2 {5^.Vn~uP{e}Vn~u5^)I{E) 

Eea{Hn)nI 

^Y. I '^(^^ + '*-(^)"') ('^x, n„„ {E) 5^) f{E) dE . (C.9) 
The eigenfunction correlator, which is defined as: 

iV„(x)i+iQo(x,y;/,s) 

' ^ ^ ' ' " (C.IO) 



V{x)^V(x)+v 

= J2 {S^^^/N^P{E}^/N^6^y'" {S^, y/NuP{E}Sy 



E£a{Hsi)nI 

can be presented in the form of ( |C.9| l, with f{E) the function which is defined in the 
neighborhood of the zeros of Xv + k^{E)^^ as 



fiE) = 



1 {S^,I^,^{E)^/Nu{Ha-E)-^Sy 



{S^,VK:P{E}Sy) 

(s^,VKP{E}VN:ils^ 



(C.ll) 



The claim follows now upon integration over v, of the expression which is obtained by 
substituting the above function in (|C.9|l. □ 



D On the Wegner estimate 

As noted in Section[2] the Wegner estimate is not being explicitly used in the Fractional 
Moment Analysis. However the statement is of intrinsic interest, and it provides also a 
useful tool for various other purposes . The basic bound has already been extended to 
multiparticle systems IICS08allKi08al . Our purpose here is to comment on the subject 
from the perspective of the approach used in Section |2] 

A local version of the bound is the following statement of finiteness of the condi- 
tional mean of the density of the spectral measure associated with the vector S-^ for an 
arbitrary configuration x G C''"'(ri), obtained by averaging over one of the potential 
variables associated with the occupied sites, V{xj; ut), j e {1, n}. In the following, 
for a self adjoint operator H we denote by Pi{H) the spectral projection associated 
with a Borel set I whose Lebesgue measure is denoted by |/|. 

Theorem D.l. Let u be a site in f2 C Z"*, x G C^"^(f2; u) a configuration with at least 
one particle at u, and 



(D.l) 
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the spectral measure of the operator Hq (lo) which is associated with the vector 
and a Borel set / C K. Then the average of fi^il', over the values of the potential 
at u satisfies: 

^ nM; ^)g{Vu)dVu ^ e (mx(/, uj) \ {v{v)}^^^) < \i\ ■ (d.2) 

Proof Bearing in mind the dependence of the Hamiltonian on Vu as expressed in ( |2.3| l, 
one may apply the averaging principle (cf. iKo86l ISW86I ICH941 ) which states that for 
any self-adjoint operator A and a bounded operator N > B ona Hilbert space: 

/ \\B^ PiiA + TN)B\\ giT)dT<\\g\\^\I\. (D.3) 

JR 

The claim, ( |D.2| i, follows from ( |D.3| l by observing that for x e C^'^^fl; u) one has: 
iV„(x) (<5x , PiiHn) <5x> = (<5x Pi{Hn) nI'^ 5^). □ 
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